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Abstract. 

We discuss the link between dark matter halos hosting the first PopIII stars formed at redshift z > 40 and the rare, massive, 
halos that are generally considered to host bright z ~ 6 quasars. We show that within the typical volume occupied by one 
bright high-z QSO the remnants of the first several thousands PopIII stars formed do not end up in the most massive halos at 
z ~ 6, but rather live in a large variety of environments. The black hole seeds planted by these very first PopIII stars can easily 
grow to M > lO^'^M© by z = 6 assuming Eddington accretion with radiative efficiency e ~ 0.1. Therefore quenching of the 
accretion is crucial to avoid an overabundance of supermassive black holes. We implement a simple feedback model for the 
growth of the seeds planted by PopIII stars and obtain a z ~ 6 BH mass function consistent with the observed QSO luminosity 
function. 
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INTRODUCTION 

Population III stars formed in the early universe at red- 
shift z > 20 with a top-heavy initial mass function (e.g. 
see fl'l) are expected to leave at the end of their lives in- 
termediate mass black remnants of the order of lOOM©. 
These seeds, formed within dark matter halos of mass 
« 10^ Mq, may be the starting points for accretion that 
will lead to supermassive black holes (Mbh > IO'^Mq), 
which are considered to power the luminosity of QSOs, 
observed in the Sloan Digital Sky Survey (SDSS) at z > 6 
when the universe was less than one billion years old 
(e.g. see (^). These bright QSOs are extremely rare ob- 
jects (one object per about 200 deg^ in SDSS, see 101), 
so we expect on average one per IGpc^ comoving. 

Within this volume the QSO may either be the descen- 
dant of the first intermediate mass black hole seed left 
from the^rif PopIII star, which would therefore give the 
most time for mass accretion, or sit at the center of the 
most massive structure at z « 6. Of course these two al- 
ternatives are in principle mutually non-exclusive, as the 
remnants of the first PopIII stars could end up in the most 
massive dark matter halos at z « 6. This possibility seems 
to be implied by a number of recent papers, where the 
progenitor halos of bright quasars are traced back in time 
and identified as the first dark matter halos formed in the 
universe (e.g. see f?], fT], fs*]). However these works ei- 
ther do not have the mass resolution to identify the dark 
matters halos hosting the first generation of PopIII stars 
or rely on multiple mesh refinements of a small region 
centered around the largest halos identified at z=0 in or- 
der to resolve scales down to lO^M©. 



To properly address the link between bright quasars 
and PopIII stars it is necessary to resolve a dynamic 
range in mass of more than 10^^: a simulation box of 1 
Gpc^ contains a mass larger than lO^^M© and within this 
box dark matter halos below lO^M^ need to be identified. 
Therefore we have adopted an original approach (see 
[6]), broadly based on the tree method by [7]. The idea is 
based on coupling a numerical simulations of structure 
formation to identify dark matter halos at z « 6 with a 
Monte Carlo method to sample subgrid fluctuations of 
the initial Gaussian random field of density fluctuations 
at the mass scale typical of halos hosting PopIII. This 
allows us to attach to every particle in the simulation, 
which has typically a mass in excess of lO'^M,::, the 
formation time of its first PopIII star progenitor. The 
details of the method are presented in detail in |6] along 
with an extensive set of tests to validate our innovative 
approach. Here we focus instead on summarizing the 
main results from our investigation. 

In particular we show that the first PopIII progenitor of 
a z = 6 bright QSO, while born at z > 40 well within the 
early era of PopIII formation, it is not the first PopIII in 
the typical Gpc^ volume occupied by the QSO, but it is 
rather preceded by about 10^ other PopIII stars. A qual- 
itative understanding can be reached from simple con- 
siderations based on the properties of Gaussian random 
fields deriving from the spectrum of primordial density 
perturbations: small mass dark matter halos are sensitive 
to higher frequency in the density fluctuations spectrum 
than their higher mass counterparts. Therefore the first 
10^ Mr.) dark matter halos formed at z > 45 in a simula- 
tion box will not in general evolve to become the first 
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FIGURE 1. From @]: Distribution of the collapse time Zni 
for the oldest PopIII progenitor (with Mj^ = lO''M0//j) of the 
most massive halo (black line) and averaged over the 10 most 
massive halos (red line) at z = 6 in the {120Mpc/h)^ box 
simulation. The blue line represents the collapse redshift of the 
first PopIII star perturbation, while the dotted green line refers 
to the collapse redshift of the 100th PopIII in the box. 



10 dark matter halos formed at z « 6. 

In terms of intermediate mass black hole growth from 
PopIII this result implies that there are a number of seeds 
formed in the early universe before the one that will be- 
come the bright z=6 QSO. All these seeds have enough 
time, if accreting at Eddington limit with accretion effi- 
ciency £ « 0.1 to become supermassive (Mbh > IO^Mq) 
by z=6. We follow their evolution and we show with a 
simple accretion model that the gas supply available for 
growth is limited for most of these seeds, so that the QSO 
luminosity function derived in our framework is consis- 
tent with the slope of the observed QSO luminosity func- 
tion. 



NUMERICAL SIMULATIONS 

We identify the largest dark matter halos at z = 6 in three 
cosmological simulations with 512^ particles and differ- 
ent box sizes: a large (edge 720 Mpc/h), a medium (edge 
512 Mpc/h) and a small (edge 60 Mpc/h) box. The sim- 
ulations have been carried out with the public version of 
the tree-PM code Gadget2 |8] and a cosmology based 
on third year WMAP data |9]: Q.a = 0.74, £2,,, = 0.26, 
H() = 70 km/s/Mpc, where O,,, is the total matter den- 
sity in units of the critical density (p^ ~ 2>HfJ{%KG)) 
with Ho being the Hubble constant (parameterized as 
//() = lOO/i km/s/Mpc) and G the Newton's gravitational 
constant.. £2a is the dark energy density. In generating 
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FIGURE 2. From |6]: Distribution for the ranking of the 
collapse epoch for the oldest PopIII halo progenitor of the 
most massive halo (black line) and averaged over the 10 most 
massive halos (red line) at z = 6 in the (llQMpc/lif' box 
simulation. 



the initial density field we use a scale invariant long-wave 
spectral index (n = 1 ) of the power spectrum of density 
fluctuations and CTg = 0.9 or CTg = 0.75 (the root mean 
squared mass fluctuation in a sphere of radius %Mpc/h 
extrapolated at z = using linear theory). As described 
in [^, the initial density field is then used as input in our 
Monte Carlo code to obtain the formation redshift of the 
first PopIII progenitor of each particle in the simulation 
box. 



RESULTS 

Under the assumption that the first PopIII stars in the uni- 
verse have formed in IQ^Mq halos cooled by molecular 
hydrogen, the progenitor of the most massive z = 6 halo 
(the one assumed to host the bright QSO) in our large 
box simulation is born at z ~ 41, while the first PopIII 
in the volume is already present at z > 49 (see Fig. [TJ. 
By the time the QSO progenitor is born, there are on av- 
erage 8000 other PopIII stars formed in the simulation 
(see Fig. |2]i, with a typical PopIII star formation rate, 
as obtained from our method, shown in Fig. [3] If PopIII 
stars reside in halos of different mass and/or if the value 
of a% is different, there will be a shift on the formation 
time of these first stars (that is on the scale of the x-axis 
in Figs. [T][3]l, but the relative ranking between the first 
PopIII star in the box and the PopIII that is the QSO pro- 
genitor remains essentially unchanged (e.g. see Fig. 3 in 

The relation between PopIII and QSOs can be easily 
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FIGURE 3. From j^: Comoving number density of virial- 
ized PopIII halos (My, = lO^Mpj//;) per unit redshift versus 
redshift for z > 29 as measured by means of our MC code. 



understood in terms of the statistics of Gaussian random 
fields (see Fig. |4]l. In fact, the most massive z = 6 halo 
in the simulation box (with mass 

originates from a « 6(7(M^/,) peak in the density pertur- 
bation field. If we know consider a random volume in 
the simulation box with mass M^/,, we find that the dis- 
tribution for the maximum peak at the mass scale of a 
PopIII halo (lO^Mo) inside this volume is in the range 
[23 : 27]a(M^/,) at the 90% confidence level. This is es- 
sentially because (i) we are at a smaller mass scale, so 
there is additional power in the density fluctuations field 
and (ii) we are considering here the first PopIII halo 
formed in the M^^/, volume, that is the maximum among 
160^ fluctuations. The probability distribution for the 
peak associated to the first PopIII progenitor of the QSO 
(green dashed line in Fig. IH is then given by combin- 
ing the probability distribution of the QSO overdensity 
with that of the first PopIII halo of a random cell. From 
Fig- El it is immediately apparent that the advantage of 
sitting at the top of the QSO overdensity is not sufficient 
for the first PopIII in this volume to beat all other PopIII 
in the box. Still one in about 10^ PopIII progenitors of 
random M^/, cells in the (720 Mpc/h)^ box is formed be- 
fore the QSO progenitor As there are « 6 x 10^ cells 
of mass Mgij in this volume, we expect from this simple 
consideration that the PopIII progenitor of the QSO will 
be formed when akeady « 6000 other PopIII stars are 
present in the box, in excellent agreement with the de- 
tailed results from the numerical simulation (see Fig.|2]l. 

This result is essentially given by the fact that bright 
QSOs at z = 6 are very rare objects. If fact, if we consider 
a smaller cosmic volume (e.g. a box of edge 60 Mpc/h) 
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FIGURE 4. From 0]: Probability distribution functions 
p{m) for the maximum m of k random Gaussian fluctuations 
representative of mass scale for halos hosting z = 6 QSO can- 
didates (red curve) and PopIII stars (black curve). To com- 
pute p{m) for QSO hosting halos we have assumed a box of 
720Mpc/h and a mass scale of Mgh = 4.3 IO'^Mq/A which 
implies k = ISO^ . This curve represents the probability distri- 
bution for a sigma peak that at z ~ 6 leads to one of the most 
massive halos in the simulation volume. The curve associated 
to first light perturbations (solid black, with Mf g = 10*M, ;i /h) is 
derived using aj^ = CT^^. — ^ and k = 160^ : it represents the 
probability distribution for the maximum of the sub-grid scale 
fluctuations at the My, scale within one cell of the 180 volume. 
The dashed green line represents the probability distribution of 
the density fluctuation associated to the first PopIII progenitor 
of a 5.7ct(M^/,) peak, m is given in units of Om^,,- 



and repeat the experiment of studying the formation time 
of the first PopIII progenitor of the most massive z = 6 
halo, we find different results (see Fig. |5]l. In this case 
the cosmic volume considered is significantly smaller, so 
sitting at the top of the most massive z = 6 halo gives 
a significant relative advantage to PopIII formation at 
z > 30: the PopIII progenitor of the most massive halo 
is typically within the first 10 to 100 PopIII in the box. 

IMBH SEEDS GROWTH 

From our investigation it is clear that, before the first 
PopIII progenitor of the most massive halo at z = 6 is 
born, several thousands of intermediate mass (niBH ~ 
IO^Mq) black hole seeds are planted by PopIII stars 
formed in a cosmic volume that will on average host a 
bright z = 6 quasar. Here we investigate with a simple 
merger tree code what is the fate of the black holes 
seeds formed up to the formation time of the quasar seed 
and what are the implications for the observed quasar 
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FIGURE 5. From d]: Like in Fig.[T] but for a smaller box 
simulation, with volume of {60Mpc/h)^ and cTg = 0.75. Given 
the limited volume available sitting on the top of the overden- 
sity that will collapse to the largest z = 6 halo gives a significant 
advantage to Poplll halos formed at z > 30, which are typically 
within the first 100 formed in the box. 

luminosity function. 

We assume Eddington accretion for the BH seeds, so 
that the evolution of the BH mass is given by: 

niBH ^mosxp[{t ~to)/t,al], (1) 

where niQ is the mass at formation time Iq and tsai is the 
Salpetertime jfioll: 

t.ai = jT^T = 4-507 • 10 Vt^ 7, (2) 

{l-e)LEdd (1-e) 

where £ is the radiative efficiency. 

If £ « 0.1 there has been enough time to build up a 
z ~ 6 supermassive black hole with mass nigH > 10^ 
starting from a PopIII remnant formed at z > 40. This 
highlights that only a minor fraction of the PopIII BH 
seeds formed before z = 40 can accrue mass with high 
efficiency, otherwise the number density of supermas- 
sive black holes at low redshift would greatly exceed 
the observational constraints. The first BH seeds in the 
box are distant from each other, so they evolve in rel- 
ative isolation, without possibly merging among them- 
selves. Therefore other mechanisms must be responsi- 
ble for quenching accretion of the first BH seeds. In- 
terestingly if we were to assume that accretion periods 
are Poisson distributed in time for each seed, we would 
not be able to explain the observed power law distribu- 
tion of BH masses at z < 6 around the high mass end. A 
Poisson distribution would in fact give too little scatter 
around the median value and a sharp (faster than expo- 



nential) decay of the displacements from the mean ac- 
creted mass. An exponential distribution of the accre- 
tion efficiency is instead required to match the observed 
BH mass function. In addition, it is necessary to assume 
that the duty cycle of the BH accretion is roughly pro- 
portional to the mass of the halo it resides in. To ex- 
plore this possibility we follow the merging history of 
PopIII halos formed at z = 40 by means of a merger- 
tree code. We implement a BH growth based on Eq. |2] 
but at each step of the tree we limit the BH mass to 
tnBH < I] nihar, where irihar is the total baryon mass of 
the halo that hosts the BH. If Tj « 6 ■ 10"^^ (like in [11]), 
then we obtain an expected mass for the BH powering 
bright z = 6 quasar of « 5 • lO^M©, which is in agreement 
with the observational constraints from SDSS quasars 
101. By fitting a power law function to the BH mass func- 
tion in the range [0.055 : 0.2] • lO'^Mir, we obtain a slope 
a « —2.6, while the slope is a w —3.7 in the mass range 
[0.2 : 1.0] • IO^^Mq, a value that is consistent within the 
1 (7 error bar with the slope of the bright end of the quasar 
luminosity function measured by [2J. 

CONCLUSION 

We study the link between the first PopIII halos collapsed 
in a simulation box and the most massive structures 
at z « 6, with the aim of establishing the relationship 
between the first intermediate mass black holes created 
in the universe and the super-massive black holes that 
power the emission of bright z — 6 quasars. We show 
that almost no correlation is present between the sites of 
formation of the first few hundred lO^'Mp//; halos and 
the most massive halos at z < 6 when the simulation box 
has an edge of several hundred Mpc. Here the PopIII 
progenitors (halos of mass Mfs ~ lO^Mg,) of massive 
halos at z < 6 formed from density peaks that are « 
1 .5a{Mfs) more common than that of the first PopIII star 
in the (720 Mpc/h)^ simulation box. These halos virialize 
around z„; ~ 40, to be compared with z„i > 48 of the^rif 
PopIII halo. 

This has important consequences. We show that, if 
bright quasars and supermassive black holes live in the 
most massive halos at z « 6, then their progenitors at 
the lO^Mp; mass scale are well within the PopIII era, 
regardless of the PopIII termination mechanism. On the 
other hand, if the ihbh / C relationship is already in place 
at z = 6, then bright quasars are not linked to the rem- 
nants of the very first intermediate mass black holes 
(IMBHs) born in the universe, as their IMBH progeni- 
tors form when already several thousands of PopIII stars 
have been created within the typical volume that hosts 
a bright z = 6 quasar The IMBH seeds planted by this 
very first PopIII stars have sufficient time to grow up to 



niBH G [0.2 : 1] • IO^'M© by z = 6 if we assume Edding- 
ton accretion with radiative efficiency £ < 0.1. Instead, 
quenching of the BH accretion is required for the seeds 
of those PopIII stars that will not end up in massive halos 
at z = 6, otherwise the number density of supermassive 
black holes would greatly exceed the observational con- 
straints. One way to obtain growth consistent with the 
observations is to limit the accreted mass at a fraction 
7] w 6 • 10^^ of the total baryon halo mass. This gives a 
slope of the BH mass function a = — 3.7 in the BH mass 
range mbh £ [0.2 : 1] • lO'^M©, which is within the la 
uncertainty of the slope of the bright end of the z = 6 
quasar luminosity function (a « —3.5) measured by |0]- 
Another important point highlighted by this study is 
that rich clusters do not preferentially host the remnants 
of the first Poplll stars. In fact the remnants of the first 
100 Pop-Ill stars in our medium sized simulation box 
(volume of {5l2Mpc/h)^) end up at z = on halos that 
have a median mass of 3 • lO'^M©//;. This suggests cau- 
tion in interpreting the results from studies that select a 
specific volume of the simulation box, like a rich cluster, 
and then progressively refine smaller and smaller regions 
with the aim of hunting for the first lights formed in the 
whole simulation (see e.g., 10. IH). Only by considering 
refinements over the complete volume of the box the rar- 
ity and the formation ranking of these progenitors can be 
correctly evaluated. 
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